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1. Introduction 

Measurement of SUSY particle ('sparticle') masses in R-parity conserving SUSY events 
at hadron colliders such as the LHC is complicated by the pair production of invisible 
Lightest Supersymmetric Particles (LSPs), one of which terminates each SUSY particle 
decay chain. The presence of these LSPs prevents direct measurement of the masses of the 
sparticles participating in the decay chains via peaks in the invariant mass distributions of 
the observable SUSY decay products. 

Several approaches to resolving this problem have been documented. These approaches 
may be grouped into three general categories. The first category consists of techniques 
which use the measurement of kinematic end-points in distributions of the invariant masses 
of combinations of visible SUSY decay products (jets, leptons etc.) from a single decay 
chain in each event. Given a long decay chain sufficient constraints may be obtained to solve 
for the individual masses ||, |||] . The second category consists of techniques which use 
the measurement of kinematic end-points in distributions of observables constructed from 
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the transverse momenta of SUSY decay products from two identical decay chains appearing 
in each event || ^, 0, These techniques can provide constraints on combinations of 
sparticle masses from very short (one-step) chains and can in principle lead to constraints 
on individual masses [0, [8|, |9], [1Q|, [llj . The third category consists of hybrid methods which 
seek to use optimally both transverse momentum and invariant mass constraints from a 
subset of events in which the same long decay chain appears in both 'legs' of each event (e.g. 



12, [L3[]). Such techniques can offer an improvement on the mass constraints obtained by 



other techniques, but may require significant integrated luminosity to generate a sufficiently 
large sample of the specific events of interest. 

This paper shall focus on the first category of mass measurement techniques - that 
involving the use of kinematic end-point constraints. The goal of the paper is to re- 
examine the constraints which may be obtained from kinematic end-points in the three- 
step sequential two-body decay chains most usually studied, with a focus on correlations 
between invariant mass observables. This will lead to a simplification of the technique 
and an improvement in the possible mass measurement precision. Similar techniques have 
been studied recently for the different case of events in which SUSY higgs particles decay 
symmetrically via two two-step sequential two-body decays [14, 15 1. Correlations between 
invariant masses were also exploited in Ref. [p!2[ . 

Section 2 commences with a review of the conventional end-point technique, while 
Section 3 re-interprets the technique from a geometric perspective. Section 4 examines 
the additional constraints which can be obtained by making use of correlations between 
invariant mass observables while Section 5 examines how these constraints may be used in 
principle to determine individual sparticle masses. Section 6 presents a case study in which 
the new techniques are used to measure sparticle masses for the SPSla benchmark model 
with fast detector simulation. Section 7 concludes and outlines some possible directions 
for future work. 



2. Kinematic constraints from three-step sequential two-body decay chains 
2.1 Definition of decay chain 

In many regions of SUSY parameter space the sparticle mass hierarchy is such that heavy 
strongly interacting sparticles produced in the initial proton-proton interaction can decay 
via a 3-step sequential two-body decay chain of the form: 

5 — > 7c — > j3bc — ► aabc, (2-1) 

where Greek letters denote SUSY states, Roman letters denote visible SM decay products 
and a is the LSP. The canonical example of such a decay chain is the decay 

<?L -»■ X2Q -> Wn<? Xih l nQ, (2-2) 

where a left-squark decays via a chain consisting of the next-to-lightest neutralino, a right- 
slepton and a lightest neutralino (LSP), in the process emitting a quark (q) and two 
opposite-sign same-flavour leptons (l n , If). The notation l n (lf) is conventionally used 
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to denote the lepton emitted nearest to (furthest from) the quark in the decay chain. In 
practice l n is indistinguishable from If which leads to ambiguities in the mass constraints 
obtained from this technique. 

In this paper we shall illustrate mass measurment techniques at both parton and 
detector- level using the SPSla benchmark SUSY model [16, 17]. The relevant mean spar- 
ticle masses for this model are m^o = 96.05 GeV, mr = 142.97 GeV, m^o = 176.82 GeV 

Xi ' Ir ' X2 

and 491.9 GeV < rrig L < 543.0 GeV. In the q~i case we assume dominant light squark 
production and a common light squark mass of 540 GeV when constructing kinematic 
bounds. The limited validity of this simplifying assumption, together with the non-zero 
widths of the squarks (especially) causes some parton-level distributions to 'leak' beyond 
the expected bounds (see e.g. Figure ^ below). For this study HERWIG 6.4 [JlS|, [l!| was 
used to generate 300k (100 fb _1 ) parton-level SUSY events filtered to require at least two 
final state leptons in each event. 



A note on notation In this paper we shall frequently refer both to global bounds on 
invariant mass distributions ('global end-points') and to bounds obtained from subsets 
of events satisfying certain selection criteria imposed on other invariant mass quantities 
('conditional end-points'). In common with previous studies we shall refer to the former 
with quantities with a single superscript 'max' or 'min'. We shall refer to the latter with 
quantities with a superscript 'max(m)' where m represents the invariant mass upon which 
the selection criteria are imposed. 



2.2 Invariant mass end-points 

The end-point technique starts by constraining combinations of masses of sparticles ap- 
pearing in the decay chain given by Eqn. ( |2.2| ) by using kinematic bounds ('end-points') 
in the distributions of the invariant masses of combinations of q, l n and If. Maxima in 
the distributions of two-body invariant masses m« are obtained when visible particles i 
and j are emitted in opposite directions in the rest frame of the \% sparticle in each chain. 
Because each decay is two-body the momenta of i and j are fixed in the rest frames of 
their respective parents. Consequently there is a one-to-one mapping between the opening 
angles of the i and j in the rest frame of the X2 anc ^ m ij- Setting this angle to its max- 
imum value of 7r radians therefore allows the position of the kinematic end-points to be 
determined analytically. The formulae are (see e.g. Ref. ]|]): 

(m 2 — m~ )(m~ — m 2 ~ ) 

Kr) = — — — — , (2.3) 

Ir 

(m~ — rn 2 ~ ) (m 2 — m~ ) 

III ~n 

,2 ^,2 \ 



lm'- — m"- )[m~ — m 2 . ) 
(m 7lT y = ij* ^ '« ^ . (2.5) 



Ir 
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For the following discussion it will also be useful to define a quantity m^ 8 * which mea- 
sures the maximum value of m q i f when the two leptons are co- linear and hence mu = 0: 



m 



max(0) \ 2 
ql f 



m 



m.r 



2 

xi 



(2.6) 



The ambiguity of assigning observed lepton momenta to l n and 1 1 conventionally pre- 
vents direct measurement of mj^ ax and m^ ax . Instead, for each event the two possible 
lepton+quark(jet) invariant masses are ordered to generate w^^U and rn 7 ^^ with analyt- 
ical end-point positions at [|| (following Ref. |||): 
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where 
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The bound provided by tti^^. arises from kinematic configurations in which m q i is max- 



imised for m q i n = m q i f = m q i^ eq y Case Al holds for the SPSla model considered through- 
out this paper. 

In addition to the two-body invariant mass end-points three-body end-points may be 
observed using m q u. The analytical formula for the maximum value of the distribution 
of m q ii is complicated by the fact that it may be generated by any one of four kinematic 
configurations, depending on the sparticle mass hierarchy. These configurations are those 
in which two of q, l n and // lie co-linear in the X2 res ^ f rame an d one contra-linear (three 
configurations) or their net momentum in the q~i rest-frame is zero (one configuration). 
Taking this ambiguity into account the maximum value of the m q u distribution is found to 
lie at § (following Ref. [|): 
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The first three cases above (B1-B3) occur if one of the three co-linear configurations pro- 
vides the maximum m q u value, with cases Bl, B2 and B3 corresponding respectively to con- 



figurations with {m q i f 
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m qin X } an< ^ { m qif = mT qt^-, m ii = m ™ aX ) m gZn = 0}- The fourth case (B4) occurs if the zero 
momentum configuration is allowed kinematically. Case Bl holds for the SPSla model 
considered throughout this paper. 

The m q n distribution also possesses a non-zero minimum when mu > 0. Conventionally 
a cut requiring mu > m™ ax /\/2 is applied, leading to a threshold in the m q u distribution 
at m^f 1 given by 1 : 




(2.10) 



Together these five mass constraints may be solved numerically [|T], |2| or analytically 
to determine the four individual sparticle masses m qL , rn^o, and m^o. In either 
case the ambiguities in the interpretation of rn^^s , and m^ ax can lead to multiple 

solutions and hence ambiguity in constraints on the underlying SUSY model. This is 
discussed extensively in e.g. Ref. m. 



3. Geometric interpretation 

The conventional approach to the derivation of kinematic end-point constraints illustrated 
above focuses on physical configurations of particle momenta which maximise or minimise 
invariant masses. The constraints are thus generated when the momentum vectors of visible 
particles are correlated in a particular manner, for instance parallel or anti-parallel. Let 
us now generalise this argument to reinterpret these correlations in a geometric context. 
Our aim will be to use these momentum correlations to identify correlations between the 
invariant masses constructed from the momenta. For convenience we shall work in the rest 
frame of the x% m the following discussion, however the calculated invariant masses are 
clearly equal to those measured in the laboratory frame. 

Our starting point is the observation that extremal values of m q u are obtained when the 
momenta of q, l n and If lie in the same plane in the rest frame of the x\- Furthermore the 
extremal values of mu, m q i n and m q i f can trivially also be obtained when this is the case. 
Neglecting the masses of q, l n and If we can write the three invariant mass combinations 
as: 

mil = 2 PinPif (1 - cos Oil) i 
m q i n = 2p qPln (l - cosflg/J, (3.2) 

m q i f = 2p q pi f (l - cos q i f ), (3.3) 



lr The value of m^JJ n is dependent on the condition imposed upon mu and therefore should strictly be 



identified in this paper as w.fJJ n , in light of the discussion in Section 2.1. We choose to retain the 
notation mJS'j" for consistency with earlier work however. 
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where pi denotes the magnitude of the three-momentum of particle i and 6ij denotes the 
opening angle between the three-momenta of particles i and j. When q, l n and If all lie in 
one plane we have the further relation: 

On + e qln + qlf = 2vr. (3.4) 

Now in the rest frame of the X2 ^ ne momentum magnitude pi n is determined from 
simple two-body kinematics to be the fixed quantity: 

2 2 
mz — m~ 

Pin = ljL - (3-5) 

The momentum magnitude p q is also fixed in this frame to be: 

2 2 
m~ — mt Q 

p q = — -. (3.6) 

X~2 



The magnitude of the momentum oilf,pi f , is however not a constant in this frame - it is 
related to the fixed momentum magn: 
Lorentz transformation dependent on 6n: 



related to the fixed momentum magnitude p| measured in the rest frame of the In by a 



2 2 
, _ m i R ~ m x? 



= 1p lf (l + pcoB6 n ), (3.7) 
where the boost factor (3 is given by: 



2 2 
mt — m~ 

0= l 2 R , (3.8) 

X2 lR 

and 7 = (1 - P 2 )- 1 / 2 . 

It is convenient now to define three dimensionless mass coordinates x, y and z: 

x= m ^ 
rriqif 



(3.9) 
(3.10) 

(3.11) 



PlnP'l f ' 

together with a dimensionless mass ratio r = m^/m-o. We can now use Eqns. ( |3.1| )- 
(|3.3| ) to express Eqn. ( |3.4[) purely in terms of invariant masses m« and momenta and then 



substitute from Eqns. (p.9j)- (|3.1l| ) to obtain a single relation between x, y and z: 

x 2 + i (1 - r 2 ) 2^ - y 2 - 2 2 ^) + (x 2 - 4)y 2 z 2 = 0. (3.12) 
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Figure 1: S 9 « surface for the SPSla benchmark SUSY model plotted in x — y — z space and 
viewed from two different perspectives. 



Eqn. (|3.12| ) defines a surface Y, q u in the 3-dimensional x — y — z space. This surface 



bounds the region within which all decay chains of the type Eqn. (^2) must lie: events in 
which the visible decay products are coplanar lie on the surface, while acoplanar events lie 
within the volume. The surface is plotted in Figure [j] for the SPSla benchmark model. A 
key feature to note is that T, q u intersects the planes with x, y or z zero along lines, since 
when any one rriij value is zero, particles i and j must be co-linear and hence rrijk must be 
uniquely determined by (and vice versa). Note in addition that if T, q u were plotted in 
m q i n — m q i f — mu space rather than in x — y — z space the radial distance from the origin 
of a point on T, q u would give m q u, and hence in these coordinates m^ ax represents the 



greatest radial distance of a point on T, q n from the origin. Eqn. ( |3.12 ) shall prove useful 



for deriving some of the kinematic bounds in two- and three-dimensions discussed below 
however it is potentially also useful in its own right - for instance for selecting events with 
the decay chain Eqn. and associating decay products with specific steps in that decay 
chain once the masses of the sparticles are known. 

4. Additional kinematic constraints 
4.1 Global end-point constraints 

Constraints on the masses of the sparticles participating in the decay chain listed in 
Eqn. ( |2.2D can be obtained by measuring the surface T, q u defined by Eqn. ( |3,12| ). The 
techniques described in Section |2.2| rely on integrating out all but one invariant mass ob- 
servable and then measuring the global end-point in the resulting mass distribution. It 
should be clear from the discussion of Section || however that there is potentially much 
more information contained in Eqn. ( |3,12 ) than can be exploited with this relatively sim- 



ple technique. In particular the shape of this surface can be measured in more detail by 
exploiting correlations between invariant mass observables. Such correlations are partially 
exploited by the definition of the ?n^ m end-point in Eqn. ( |2.10| ), in which the integration 
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Figure 2: Parton-level two-dimensional invariant mass-squared distribution for decay chain 



Eqn. (2.2) for the SPSla benchmark model showing m 2 qll 



plotted against m^. See text 



for explanation of bounds. 



takes place only over a limited range in mu, however the possibilities presented by this 
technique are more numerous, as shall be illustrated in the following discussion. 

4.2 Conditional end-point constraints with two-dimensional correlations 

The first class of additional constraints can be obtained by integrating over only one, rather 



than two, of the three degrees-of-freedom used to define the surface H q u in Eqn. ( 3.12| ). 
The most familiar example of such a technique is that used to obtain the m^f 1 end- 
point discussed above. The correlation between m q u and mu can be observed by plotting 
g K — rnfi against m u , shown at parton-level in Figure |2| for the SPSla benchmark model. 



m 



The vertical kinematic bound in this figure is provided by m™*™ (right-hand vertical line), 
while the minimum of the one-dimensional distribution obtained by integrating along the 
mfi axis to the right of the left-hand vertical line measures (m^f 1 ) 2 — (m™ ax ) 2 /2. 

The curved bounds in Figure Q may be obtained from the general formula (see e.g. 
Ref. |20| 1) for the end-points in the invariant mass distribution of two heavy visible decay 
products a and b of a two-step sequential two-body decay chain 7 — > fib —> aab: 
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bound (m a ,m,f,) \ 2 



ah 



) =m a +m b + 
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(4.2) 
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When a = 11 and b = q and mj, = this equation may be simplified to give: 



( 



m 



qll 



m U + 



m 



A2 



«) 



2m 2 



m 2 ~ 



m 2 ~ + 



m\ ± ^Xim^m^mu)^ , (4.3) 



where the positive root corresponds to the usual expression for the 'm™^ x ' end-point with 
X = 11 HI]. The two roots of this equation correspond to the upper and lower curved 
bounds in Figure ||. Setting m a = with a = If and b = ql n in Eqn. ( |4.1| ) also allows the 
equivalent bounds in the m 2 u 



m 



qln 



versus rn q i n distribution to be obtained (not shown). 



Plots similar to Figure |] can also be constructed from the distributions of m qll 



m 



as a function of m 2 , and m' 



gin 



2 , as a function of m 2 , . In the latter case however 



2 

l qll ~ m qlf ao a ^^"^ "i. ,,o qlf . 

the kinematic bounds are particularly difficult to describe analytically. An alternative 
means of illustrating two-dimensional correlations in which the kinematic bounds may be 
described more easily is provided by the distributions of m?- as functions of m 2 k where i, j 
and k are drawn from q, l n and If. The three possible distributions are shown in Figure || 
at parton-level for the SPSla model. Squared masses are plotted in this and subsequent 
figures in this section (as in Figure |2|) as this generates kinematic bounds which take on 
particularly simple forms. 

In Figure || the bounds on m 2 ,^, m qi n an d mfi obtained from the conventional 'global 
end-point' analysis together with ( m ^ X ^) 2 are illustrated by vertical and horizontal lines 
(labelled). In the top-right and bottom figures the upper bounds on m 2 , . lie respectively 



at: 
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ql f 



\^ / max(0)\2 , 
) = K m al f ) + ( 
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(4.4) 



and 
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max(m gi ,J^2 
ql f 



) =WD - Kr) -( 
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max(0) N 1 2 
ql f 
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m qlr, 



(4.5) 



Of course in practice we can not distinguish m q i f from m q i n on an event-by-event basis. 
For this reason we plot in Figure the equivalent distributions to Figure 



m 



rather than m 2 . 



and m 2 qlf . 



for m2 ql(lo) and 

As expected the distributions resemble combinations of 



<?Z(hi) ""ql n 

the rn q i n an d m qi f distributions, with an additional bound where m q i n = m q i f = m q ^ eq y 
This bound is trivial in Figure |3j(bottom) while in the top-left and top-right figures the 
equation of the bound can be obtained by using Eqn. ( 3.12 ) and Eqns. (|3.g|)-( |3~Tl] ) and 



solving the resulting quadratic equation for Tn 2 ^^ as a function of ruff. 



H + A 2 ) 



2\ 4 
m ql(eq) 



2m llPq (2pi n + A)m ql{eq) + m u p q = 0, 



(4.6) 



where 



A^rp' lf -p ln + (l-r 2 )^ 



(4.7) 



The two roots of this equation give the upper and lower curved bounds in these figures. 
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Figure 3: Parton-level two-dimensional invariant mass-squared distributions for decay chain 



Eqn. (2.2) for the SPSla benchmark model. The top-left figure shows the distribution of m ql as 



a function of mfi , the top-right figure the distribution of rn qlf as a function of mf[ and the bottom 
figure the distribution of ra 2 q i f as a function of m qln . See text for explanation of bounds. 
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Figure 4: Observable two-dimensional parton-level invariant mass-squared distributions for decay 
chain Eqn. (2.2) for the SPSla benchmark model. The top-left figure shows the distribution of 
m qi{io) as a function of mfi, the top-right figure the distribution of rn? ql ^ as a function of mfi and 
the bottom figure the distribution of m 2 ql ^ as a function of m q u\ y See text for explanation of 
bounds. 
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Figure 5: Two-dimensional distributions of two-particle invariant masses for the SPSla model, 
with the third invariant mass required to be less than 50 GeV (top figures) or less than 10 GeV 
(bottom figure). The top- left figure shows vrr qln versus m^, the top-right figure ra q i versus mj^ 
and the bottom figure versus . See text for explanation of bounds and loci. 



4.3 Conditional end-point constraints with three-dimensional correlations 



A second class of constraints can be obtained from Eqn ( p. 12 ) by measuring invariant mass 



correlations without integration over additional mass distributions. In principle this could 
be used to measure T, q n directly. Here we provide one example of such a technique in 
which we study the form of the one-dimensional lines where S„^ intersects the planes with 
respectively x, y and z (or equivalently m q i n , m q i f and ma) equal to zero. The resulting 
mass-squared distributions are analogous to those shown in Figures | and | but with the 
additional requirement that the third mass(-squared) value not plotted in each figure is 
small. 

We proceed by requiring m q i n , m q i f or m« to be less than some value (here set to 50 
GeV for m q i quantities and 10 GeV for mu) 2 sufficiently large to provide useful statistics but 



2 With these selection requirements the mean quark-lepton(lepton-lepton) separation is AR = 0.6(0.2). 
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Figure 6: Observable distributions constructed from the quantities plotted in Figure |5| for the 
SPSla benchmark model. The lower of the two diagonal lines in the right-hand figure represents 
the trivial bound where m q i^ — m q i(\ y The remaining bounds and loci are as for Figure [5] (see 
text). 



sufficiently small to approximately maintain the one-dimensional form of the intersection. 
We then plot the two-dimensional distribution of the two unconstrained squared masses in 
each case. These distributions can be seen in Figure [51 In the top-left and top-right figures 



the vertical lines represent the bound from (m™ ax ) , while the horizontal lines represent 

"qln ) > V'-qlf J l"V/ 



the bounds from (w-^ ax ) 2 , ( m q c i a ' x ) 2 an d i 112 ^^^) 2 • ^ n ^ ne bottom figure the vertical and 



horizontal lines represent the bounds from (w^ x ) 2 and (w^J*^) 2 respectively. The loci 
of the distributions in the figures are given in terms of dimensionless mass coordinates by: 

x 2 (y = o) = — rrr 2^ 2 » ( 4 - 8 ) 

r + 4 (1 — r z ) z £ 

y \ x = 0) = z 2 , (4.9) 

y 2 (z = 0) = x 2 r, (4.10) 

which gives the following relations: 




< = ml [ -^-T^I \ , (4-11) 



<l f = < I _2 J 2 I , (4-12) 
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™% S = m 2 qln \ J« I • (4-13) 



If these object separations were too small in practice to provide useful lepton reconstruction efficiency they 
could be increased at the expense of mass resolution (see text). 
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Figure 7: Distribution of the ratio ^gi(hi) / m g/(io) f° r m u < 10 GeV for the SPSla benchmark 
SUSY model. The distribution possesses a trivial bound at rriqi^ / m q i(\ ) — 1 while the prominent 
peak is indicative of sequential two-body lepton producing decays. 



The distributions in Figure || are not observable directly due to the l n /lf ambiguity. 
Instead the three distributions of Figure || are replaced by two distributions (with m q i(\ ^ 
required to be small), plotted in Figure ||. If the loci in this figure can be identified 
in the presence of realistic detector-level smearing then they could be used to measure 
sparticle masses via Eqns. (4.11)~( |4.13| ). In addition however Figure ||(right) can be used 
to determine whether the decay chain in selected events consists of two sequential two-body 
lepton-producing decays (as in Eqn. ( |2.2[ )) or a single three- (or more) body decay (e.g. 
X2 - * l + l~Xi)- I n the latter case mjj can acquire a small value through one of the two 
leptons acquiring a very small momentum in the X2 res t frame, while in the former case 
this can only occur due to event topology, i.e. due to the two leptons being emitted co- 
linearly in this frame. This conclusion may also be reached by observing that T, q u defined 
in Eqn. ( 3,12| ) does not intersect any of the x, y or z coordinate axes and instead intersects 
the x = 0, y = and z = planes in one-dimensional lines generating the loci shown in 
Figure |H[ The consequence of this observation is that in the three-body case m^^i) and 
m qi(io) should be rather less correlated than in the sequential two-body case, where m q in^\ 
is linearly related to m q m \ as shown in Figure |6|(right). A convenient means of observing 
this correlation is to plot the one dimensional distribution of the ratio m q i(hi) / m qi(io) shown 
in Figure [/], where the presence of the prominent peak provides evidence for the presence 
of sequential two-body lepton producing decays. This provides a useful alternative to 
conventional techniques such as studying the shape of the mu distribution or measuring 
the ratio (which equals y/2 for three-body lepton-producing decays pif] )- 
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The kinematic bounds for three-body lepton-producing decay chains incorporating single 
two-body quark-producing decays are listed for completeness in Appendix A. 



5. Mass reconstruction 

Having identified some correlations between invariant mass observables we shall now discuss 
how they might be used in principle to constrain individual sparticle masses. 

We start by observing that measurements of the global end-points together with mea- 
surements of the straight-line upper bounds in the two-dimensional distributions shown in 
Figure || can in principle give direct access to the four 'core' two-body observables m™ ax , 
m^ x , tu^ x and ti^J* . In order to measure these quantities however we must first re- 
solve or evade the ambiguities in the interpretation of the m q \ and m q u end-points listed in 
Eqns. ( |2.7D and fl2.9| ). This can be accomplished by examining the form of the kinematic 
bounds in Figures |and|. The upper m 2 qln bound (top-left panel of Figure |||) is constant as 
a function of m u , while the equivalent bound (top-right panel of Figure |3j) is a linear 
function of mf, with non-zero gradient (in the absence of sparticle mass degeneracies). Con- 



sequently if the conditional upper bound on m^( m ) changes with win then m'7/C: ~ '"'! 



.max m max 

V(hi) - m qlf ' 

corresponding to cases Al and A2 in Eqn. fl2.7[ ). This is the situation for SPSla. Con- 
versely if the upper bound on m q i^ remains constant with mu then m^r^^ = m ^ x > which 
corresponds to case A3 in Eqn. ( |2.7D , and hence also tn^a^ = rn^^y 

In order to discriminate between cases Al and A2 we could in principle examine the 
behaviour of the conditional upper bound on m q i(\ ^ as a function of mu- If the bound 
is constant in the vicinity of Tn^^s then tti^q^n = "i^ ax (case Al, realised at SPSla), 
while if it changes as a function of mu then n^^S = m^?^ (case A2). Equivalently, 
if the distribution of m g u m ) values as a function of m q i^ (shown for squared masses in 
Figure ^(bottom)) is bounded in m q i(i ^ by a vertical line then rn™^^ = m™j^ (case Al) but 
if it is bounded by a point then Tn^y\ = tit^^. (case A2). Nevertheless the (squared mass) 
distributions shown in Figure |i] for SPSla, for which m™^ is very close to Tn^?^., illustrate 
the difficulty of this strategy for resolving the tti^^, ambiguity. It is far from clear in these 
parton-level distributions whether these conditions are satisfied and at detector-level the 
difficulties will be still greater. In this case however we can potentially also make use of 
the fact that the upper bound on rn? qln shown in Figure ^(top-left) is also visible in the 
top-right and bottom panels of Figure [| as a discontinuity in the ^^i) distribution below 
its upper bound. Consequently observation of a 'hidden' bound in the m^(hi) distribution, 
with a position independent of any conditions imposed upon mu would enable w-"J ax to be 
measured directly, evading the A1/A2 ambiguity. 

In cases Al (e.g. SPSla) or A2 the conditional upper bound on m q i^ as a function 

of mu allows w-^ ax and fn^^ to be determined with Eqns. (|4,4|). Together with the 
measurement of ?n<^ ax discussed above and the measurement of the m™ ax end-point all 
four core two-body observables can hence be accessed. In case A3 it is the upper bound 
on m q i f which is hidden but nevertheless it may still be observable below the upper bound 
on m qi ^ provided here by the upper bound on m q i n . If this is the case then all four 
core two-body observables can again be determined. There are certainly special cases 
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where the sparticle mass values conspire to prevent ambiguity resolution on the basis of 
these arguments, however in those cases the use of multiple redundant measurements, for 
instance bounds on m q u, should help to resolve these ambiguities (see e.g. Ref. [f|). The 
mu dependence of the conditional upper or lower bounds on m q in^\ an d Tn q m \ (shown for 
squared masses in Figure ||) could also be exploited. 

An additional input to the mass reconstruction can be provided by the global upper 
bound on m q u provided by m^J ax . This bound also suffers from ambiguities, as listed in 



Eqn. (2J}). Having resolved the ambiguities in m q i however these ambiguities are consider- 



ably easier to address. The three cases corresponding to co-linear kinematic configurations 



(B1-B3) possess (m^ ax ) values of respectively {m q ^ f ) + ( m ^ x ) > ( m ™ ax ) +( 
and (m^ ax ) 2 + (m™ ax ) 2 . Consequently the position of the value of m 



max 
qll 



11 J ^V"V 

is already uniquely 

2 



specified by independent measurements of the four core two-body observables. If (wi^f x ) 
does not equal any of these three values within errors then case B4 is likely to provide 
the bound and hence an additional kinematic constraint. It should be noted that further 
mass constraints can in principle be obtained using Eqn. ( |2.9D by measuring the depen- 
dence of the conditional upper bound on m q u as a function of mu together with the similar 
dependence of the conditional lower bound on m q u, from which the global m^f 1 end-point 
derives. 

If we can unambiguously measure at least the four core two-body observables then we 
can obtain the individual sparticle masses analytically from: 
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(5.1) 

(5.2) 
(5.3) 
(5.4) 



Additional constraints, for instance provided by "i^ ax , f^/ 111 or measurements of the con- 
ditional upper or lower bounds on m^hi) and m q uy\ can be included in several equivalent 



ways. One approach involves replacing the analytical formulae of Eqns. (5.1)— (5.4) with 
a numerical fit to the sparticle masses incorporating all constraints. This approach is dis- 
cussed further in Section |(| An alternative procedure involves re-casting the constraint 
equations in terms of the four core two-body observables and using a numerical fit to de- 
termine these quantities, retaining the analytical formulae for determining the sparticle 



masses in a second step. In this approach Eqns. ( |4.3| ) and ( [4.6[ ) can be re-cast in terms of 
'"' xt0 ^ and m[^ ax with Eqns. (5.1)-( pTi| ) together with the equations relating 
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these quantities to the momenta of the jets and leptons considered in Section ^ 

ty-y% nicLX/v-v-l m£*x 

i it'll qln 

Pin = 2m max" > ( 5 ' 5 ) 
ql f 

ql ql 

Pq = , (5-6) 

^ = • ( 5 - 7 ) 

6. Detector-level study 

6.1 Introduction 

Having outlined some possible sparticle mass measurement strategies making use of the 
additional information provided by invariant mass correlations we now illustrate how such 
strategies may be used in practice with detector-level events. Although a very wide range 
of new variables are potentially now available we shall concentrate on a small subset repre- 
senting a limited evolution of the conventional global end-point technique. This will enable 
us to highlight the main benefit of the new technique, namely the amiguity resolution 
discussed in Section [|, without compromising event statistics near the end-points. 

We shall focus purely on m q u and m q i end-points as functions of mu in order to 
minimise the possibility of under-estimating statistical uncertainties through the presence 
of correlations between end-point observables. Such correlations could arise if the same 
events appear at the end-points of two or more invariant mass distributions. Estimation 
of these correlations is beyond the scope of the current paper and so we choose instead 
to minimise their effects by judicious selection of minimally correlated observables. With 
more work to understand such correlations further improvements in mass measurement 
precision could potentially be obtained. 

It should be noted that throughout this section we consider end-points in distributions 
of invariant masses rather than squared invariant masses. The use of the former here 
enables comparison of results with those of previous studies, and is consistent with the 
discussion of Section ||. 

6.2 Event simulation, selection and reconstruction 

The 100 fb _1 equivalent sample of SPSla events described in Section ^]l] was passed through 
a fast simulation of a generic LHC detector [22]. In addition a 100 fb _1 equivalent sample 



of tt background events was generated with HERWIG 6.4 [18], passed through the same 
fast detector simulation, filtered to require at least two leptons (e or fj,) and added to the 
signal sample, tt events are expected to form the dominant SM background to the analysis 
described here, however the event selection and background subtraction described below 
are found to reduce this contribution to negligible levels (typically < 10 events for 100 fb~ x 
of data). 



Events were selected with the same requirements as were used in Ref. [17] to aid 
comparison: 
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• At least four jets with the hardest three satisfying: p-r(ji) > 150 GeV, pt{J2) > 100 
GeV, PT (h) > 50 GeV 

. M eS = E™ ss + PT (ji) +Pt{32) +Pr(h) +Pt(h) > 600 GeV 

• £f? iss > max(100GeV,0.2M eff ) 

• Exactly two isolated opposite-sign same-flavour (OSSF) leptons (e/^) satisfying: 
p T {h) > 20 GeV, pr{h) > 10 GeV 

Following event selection the invariant mass of the OSSF lepton pair (mu) was calculated 
for each event. The mu distribution of events containing opposite-sign opposite-flavour 
(OSOF) lepton pairs was subtracted from that of OSSF events in order to remove SM and 
SUSY background events containing uncorrelated leptons. The resulting mu distribution 
is shown in Figure |8|(top). This OSOF subtraction procedure was also applied when gen- 
erating all of the invariant mass distributions shown below. These distributions also all 
require that mu < 80 GeV (3 GeV above the observed value of m™ ax ). 

The next step in the analysis involves indentifying the hard jet generated by the as- 



sumed decay chain Eqn. (2.2). Throughout this analysis this jet was assumed to be one 
of the two leading jets in the event. When calculating the maxima of invariant mass dis- 
tributions this jet was further assumed conservatively to be the one which minimises m q u, 
while when calculating thresholds and invariant mass ratios the jet which maximises m q u 
was used. This latter m q u value is referred to as m q un^\ below. When calculating m q i 
end-points only those events in which the two possible m q u values lie either side of the 
observed m^f x end-point were used. When calculating distributions of ratios of invariant 
masses this requirement was strengthened to require that both m q u values lie below the 
observed m^ ax end-point. 

6.3 Global end-point analysis 

The integrated m q u, m q un£), fn q u\ ) an d m q i(hi) distributions obtained from this analysis 
are shown in the bottom four panels of Figure M and agree well with the results of earlier 



work [17]. These distributions were fitted with linear functions to determine the posi- 
tions of the global end-points while the mu distribution (Figure |8](top)) was fitted with a 
'triangular' distribution smeared with a gaussian resolution function of variable width, as 
described in e.g. Ref. Q. A more detailed analysis might involve fitting with the analytical 
functions described in Ref. Q or simulated detector-level template distributions, however 
such developments are beyond the scope of this paper. 

The mean end-point positions and their associated statistical uncertainties obtained 
from the fits are shown in Table |l]. The agreement between the nominal end-point positions 
expected from the masses used in the HERWIG generator and the fitted means is reasonable 
in all cases. Also shown in Table [l] are the jet/lepton energy scale (JES and LES) system- 
atic uncertainties, assumed to be 1% and 0.1% respectively. The scale uncertainties in the 
m q i and m q u end-points are dominated by the JES uncertainty leading to a 100% corre- 



lation between the systematic uncertainties in these observables. As discussed in Ref. [17] 



- 18 - 



a 

o 
o 

> 

o 



800 
600 
400 
200 




v> tt 



tttt 



t + ++ 



10 20 30 40 50 60 70 80 90 100 

m u (GeV) 



o 
o 



> 

o 



o 
o 




£ 300 



200 



100 



400 
300 
200 
100 



m || >m 1| max/V2 



400 600 

m qll (GeV) 




400 



300 



200 



400 600 

m qll(hi) ( GeV > 



200 400 600 

™ q i,.„)(GeV) 




400 600 

m ql(hi) (GeV) 



Figure 8: Dctcctor-level invariant mass distributions for the SPSla benchmark SUSY model used 
in the conventional integrated end-point analysis. 

the systematic uncertainties arising from variation of input parameters in these naive fits 
are potentially significant (<10 GeV) however for the purposes of this analysis we shall 
henceforth neglect these systematics. 
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Observable 


Input 


Mean 


Statistical error 


Energy scale error 


max 

"hi 


77.053 


77.006 


0.057 


0.077 


m max 


428.4 


425.3 


1.8 


4.3 


min 
" l qll 


201.8 


200.8 


3.1 


2.0 


max 
m ql(hi) 


377.9 


378.0 


1.6 


3.8 


max 
m ql(lo) 


300.2 


301.3 


1.0 


3.0 



Table 1: Summary of measurements of end-points in integrated invariant mass distributions for 
the SPSla benchmark SUSY model, states. Column 2 lists the end-point positions expected from 
the masses used in the HERWIG generator while Columns 3 and 4 provide the fitted end-point 
positions. Column 5 provides the systematic uncertainty obtained from jet and lepton energy scale 
uncertainties of 1% and 0.1% respectively. All masses are in GeV. 



In order to determine the precision with which individual sparticle masses can be 
measured with the global end-point constraints listed above, 10000 LHC experiments were 
simulated with a simple toy Monte Carlo code. For each MC experiment values for each 
of the end-point positions were drawn from gaussian distributions centred on the nominal 
values with widths determined by the uncorrelated statistical and correlated systematic 
scale uncertainties listed in Table |l|. The end-point positions were then fitted with the 
formulae listed in Section U (Eqns. (fl|), (^7|) ; ( ^9|) and (ETIc| )) using MINUIT ||. The 
X 2 minimisation function which was used was that described in Refs. [||, |j: 

x 2 = j E exp _ E th (m)] T W[E cxp _ E th (m)], (6.1) 

where E exp and E th (m) are vectors of respectively observed and predicted end-point posi- 
tions (the latter functions of the vector of sparticle masses m), and the weight matrix W 
is the inverse of the correlation matrix of the end-point observables. This form of the min- 
imisation function takes into account, through the use of W, the correlations between m q u 
and m q i end-point observables generated by the correlated JES systematic uncertainties. 

As discussed in Ref . j|] , excluding the m^f 1 constraint generates a two-fold ambiguity 
in the allowed set of mass values for the SPSla model, corresponding to cases Al and A2 
in Eqn. (^?j). To properly account for this ambiguity we performed the fit twice for each 
MC experiment, with the initial values of the four masses set to the analytical solutions 
given by these two cases, listed in Ref. |J. The solution with the minimum x 2 value was 
then selected. 

The distributions of fitted mass values obtained from the 10000 toy MC experiments 
are shown in Figure |9], which also provides the means and RMS values of the distributions. 
The latter provide estimates of the sparticle mass measurement precisions. Also shown in 
the Figure are the distributions of mass solutions obtained by starting the fits from each 
of the two sets of analytical formulae corresponding to cases Al and A2. The incorrect 
solution corresponding to case A2 contributes 32% of solutions and generates mass values 
~ 10 GeV below the nominal values. It can clearly be seen that these incorrect solutions 
significantly decrease the precision with which the masses can be measured. 



- 20 - 



2 



W 

u 



e 
Z 



1200 
1000 
800 
600 
400 
200 




60 80 



100 120 140 

m( Xl ) (GeV) 



100 120 140 160 180 200 

m(l R ) (GeV) 




140 160 180 200 220 

m(x 2 ) (GeV) 



500 520 540 



560 580 600 

m(q L ) (GeV) 



Figure 9: Fitted sparticle masses for the SPSla benchmark SUSY model obtained from the 
conventional global end-point analysis. In each case the two smaller histograms underlying the 
large (sum) histogram represent the distributions obtained by using the analytical mass formulae 
for cases A2 (left - incorrect assumption for SPSla) and Al (right - correct assumption for SPSla). 
The latter distributions possess widths ranging from 6.6 GeV (m^o) to 11.3 GeV (rriq L ). 



6.4 Conditional end-point analysis 

When exploiting invariant mass correlations with conditional end-points we start by con- 
firming that the selected events do indeed contain sequential two-body lepton-producing 
decays. In the global end-point analysis this can be accomplished by studying the shape of 
the mu distribution (Figure ||(top)), which departs from the canonical 'smeared triangle' 
when the dilepton pair is produced off mass-shell (see e.g. Ref. |2l|). This shape analysis 
potentially requires large event statistics and a good understanding of detector performance 
to conclusively exclude the three-body hypothesis however, since the latter is effectively 
the off mass-shell limit of the sequential two-body case. With invariant mass correlations 
however we can also plot the reconstructed ratio m q i(hi) / m qi(io) Ior small mu discussed in 
Section 4.3. This is plotted in Figure [H] at detector-level for mu < 10 GeV. Comparison 
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Figure 10: Distribution of the detector-level ratio m g ;(i 1 i)/ TO gi(io) f° r m u < 10 GeV for the SPSla 
benchmark SUSY model. 



with Figure [7| shows that the prominent peak associated with sequential two-body lepton- 
producing decays has survived detector-level smearing and the application of selection cuts 
and hence this distribution provides clear evidence for such decay chains. 

The next step in the analysis involves identifying whether m^^^ measures m^ x or 
m^ x , i.e. whether cases A1/A2 or A3 in Eqn. ( |2.7| ) are correct. To accomplish this we 
construct two m^Oji) distributions (Figure |ll|) - one for events with low mjj (< 30 GeV) 
and one for the remaining events (30 GeV < mu < m™ ax ). Fitting the end-points ^^n" 1 "^ 
of these distributions with simple linear functions we find that their positions differ signif- 
icantly (see Table Q). This confirms that the value of w^mj\ used in the global end-point 

analysis measures m^ x , as does the larger of the two m ^rbd) end-points measured here. 
The smaller of the two ^^ry) end-points can then be used with Eqn. ( |4.4| ) to determine 
111,1 • sll) ) as described in Section |5[ 

Having excluded case A3 from Eqn. ( |2.7|) we now aim to measure m^ x unambiguously, 
which will help us to resolve the A1/A2 ambiguity in rn^^y To do this we construct the 
distribution of m q i, where each selected event contributes two entries (m^uy) and m q i^). 
We further require that mu be large (70 GeV < mu < m™ 8 *) in order to maximise the 
separation between the conditional upper bound on m^hi) provided by m^**" 1 "-*, at the 
lower end of the mu range, and the 'hidden' end-point measuring m™^, whose position is 

independent of mu- The value of m^ ax ^ m at the lower end of the mu range determines 
the value of m q i at which the distribution begins to fall off ~ linearly towards m^ ax . 
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Figure 11: Detector-level m^^i) distributions for the SPSla benchmark SUSY model for mu < 
30 GeV (left) and m u > 30 GeV (right). 
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Figure 12: Detector-level m q i distribution for the SPSla benchmark SUSY model for mu > 
70 GeV (left) and showing end-point fit (right). The same distribution is plotted in the left-hand 
figure as the right-hand figure but with twice the bin-width to improve the visibility of the m q i n 
end-point. 



Consequently by maximising the separation of this m ™^ x ( mn ) value from the end- 
point we maximise the visibility of the latter. The resulting m q i distribution is shown in 

380 

end-point at 



Figure 12 and displays both the linear end-point expected from m"J ax ^ m "^ at 340 



GeV and the hidden m" 



300 GeV. The lack of dependence of the m 



max 
qln 
max 



end-point on mu causes it to have the same characteristic 'triangular' shape as the m u 
end-point at parton-level and consequently we fit this distribution with the same 'smeared 
triangle' used to determine the m™ ax end-point position. The fit results are shown in 
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Observable 


m H (low) 


m«(bigh) 


Input 


Mean 


Statistical error 


Energy scale error 


max 
"Hi 






77.053 


77.006 


0.057 


0.077 


max(mn) 
'"qll 


0.0 


70.0 
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1.4 
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IllelXl 'ATI// / 

m q U 
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400.0 


3.7 


4.0 


™min 
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™ max / /o 
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ZUl.O 


zUU.o 


Q 1 
6.1 


z.U 


max(mn) 
m ql(hi) 


0.0 


30.0 


317.6 


322.1 


3.2 


3.2 


max(mn) 

max 
" l ql n 


30.0 
70.0 


max 

"hi 

lyy-i max 

"hi 


377.9 
300.2 


379.6 
295.4 


1.8 
2.7 


3.8 
3.0 


max(mii) 

"vo°) 


0.0 


70.0 


300.2 


300.5 


1.2 


3.0 


max(mii) 

m ?;(io) 


70.0 


^y-y-j max 

m ll 


268.3 


268.1 


2.6 


2.7 



Table 2: Summary of measurements of end-points in conditional invariant mass distributions for 
the SPSla benchmark SUSY model. Columns 2 and 3 provide the mu ranges over which the 
invariant mass distributions are integrated to give the conditional end-points listed in Column 1. 
Column 4 lists the end-point positions expected from the masses used in the HERWIG generator while 
Columns 5 and 6 provide the fitted end-point positions and statistical uncertainaties. Column 7 
provides the systematic uncertainty obtained from jet and lepton energy scale uncertainties of 1% 
and 0.1% respectively. All masses are in GeV. 



Table ||. The visibility of the hidden m™j^ end-point may be less clear if a more realistic 
detector simulation is used, however more detailed studies of the m q i distribution over the 
full mu range may corroborate an ambiguous observation in this case. Henceforth we shall 
therefore assume that an unambiguous measurement can be obtained. 

With the extra information provided by the ni^ x end-point we should potentially be 
able to resolve the A1/A2 ambiguity in m^^y Rather than simply measuring the value 
of the global m-g^io) bound however we proceed by measuring the value of the conditional 
bound w-^o)™"^ hi two mu bins: < mu < 70 GeV and 70 GeV < mu < m™ ax . The first 
bin contains the global end-point w-^^ (which for SPSla measures m™f^ via case Al) 
while the second provides a measure of the shape of the upper bound on m q u\ Q ^ determined 
by Eqn. (|4,6|). The m q i(\ -j distributions in these bins are shown in Figure 13 and the 
results of linear end-point fits are listed in Table Of course we do not know a priori how 
to interpret these end-points, however by using these measurements as constraints in the 
subsequent global mass fit the correct interpretation can be determined with the help of 
the other end-point measurements, particularly m™^. 

We also measure the conditional bound "rn^^™ 11 ^ in the same two mu bins used above 

to measure m ™i^ an ' '• The two m q u distributions are shown in Figure 14 while the results 

of linear end-point fits are listed in Table |2[ The fn^^^ 11 ^ end-point observed in the first 
bin lies at the position of the global m q u end-point w.^f x while the end-point in the second 

bin provides a measure of "rn^f^" 111 ^ determined by Eqn. (|4[ 

In addition to the above end-point constraints we make use of the m™ ax and m^lf 1 
constraints from the global end-point analysis. When added to the end-point constraints 
discussed above this gives nine end-point constraints, listed in Table |2[ In most cases there 
is reasonable agreement between the fitted means and the nominal end-point positions. 
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Figure 13: Detector-level m 9 ;n ) distributions for the SPSla benchmark SUSY model for mu < 
70 GeV (left) and m u > 70 GeV (right). 
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Figure 14: Detector-level m q u distributions for the SPSla benchmark SUSY model for mu < 70 
GeV (left) and m u > 70 GeV (right). 



The fitted position of the first n<^f x ^ m "^ end-point departs from the nominal position by 
4.6a indicating that more work is needed to define the fitting function in this case, however 
this will not be considered further here. 

As in the global end-point analysis the precision with which SUSY particle masses 
can be measured was determined by generating 10000 toy MC experiments in which the 
end-point positions were smeared from their nominal values with gaussians with widths 
determined by the appropriate uncorrelated statistical and correlated systematic scale un- 
certainties. A global fit to the end-point constraints was then performed with the formulae 
from Sections 2.2 and |4.2| as described in Section |6.3|. The initial values of the sparticle 



- 25 - 



% 1400 



a 

1) 



W 



e 
Z 



1200 
1000 
800 
600 
400 
200 



1600 




100 120 140 

m( Xl ) (GeV) 



100 120 140 160 180 200 

m(l R ) (GeV) 



900 




140 160 180 200 220 

m(x 2 ) (GeV) 



500 520 540 



560 580 600 

m(q L ) (GeV) 



Figure 15: Fitted sparticle masses for the SPSla benchmark SUSY model obtained from the 
conditional end-point analysis using invariant mass correlations. 

masses used in the fits were determined by calculating m^ x (g iven by the larger n^qiQ^ u ^ 
end-point), Jn^** (obtained from the two rr ^^^ ni1 ^ end-points), rn™j^ and m™ ax and 



then using Eqns. (|5.1[) - (|5.4j) . 

The distributions of sparticle masses obtained from the 10000 MC experiments are 
listed in Figure [l5|. The means and RMS values of these distributions are compared with 
the equivalent quantities from the global end-point analysis in Table [3|. The distributions 
obtained from the conditional end-point analysis are ~ 20%-30% narrower, and also less 
skewed, than those obtained using global end-points, primarily due to resolution of the 
A1/A2 ambiguity but also due to the use of additional end-point constraints. It should be 
noted however that as in the case of the global end-point analysis the conditional end-point 
analysis does not establish the identities of the intermediate states, merely their masses. 

7. Conclusions 

This paper has shown that additional constraints on sparticle masses at the LHC can be 
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State 


Input 


Global end-point fit 


Conditional end-point fit 


Mean 


Error 


Mean 


Error 


X? 


96.05 


94.8 


8.2 


95.8 


6.1 


Ir 


142.97 


141.4 


8.8 


142.7 


6.1 


xl 


176.81 


175.4 


8.3 


176.5 


6.0 


QL 


540.0 


538.5 


12.8 


539.3 


10.2 



Table 3: Summary of sparticle mass measurement precisions for SPSla states. Column 2 lists 
input masses used in the toy MC simulation, Columns 3 and 4 the fitted masses and uncertainties 
obtained from the conventional global end-point analysis and Columns 5 and 6 the equivalent values 
obtained from the conditional end-point fit. All masses are in GeV. 

obtained by exploiting correlations between invariant mass end-points obtained from the 
sequential two-body decay chains used in the conventional global end-point analysis. The 
same techniques can also be used to confirm that selected events contain sequential two- 
body lepton-producing decays rather than three-body decays. These additional constraints 
can be used to improve the precision with which sparticle masses can be measured, although 
care must be taken to acccount for any correlations between end-point positions. 

The results of this paper indicate a number of potential avenues for future work. Most 
importantly a detailed detector-level Monte Carlo study of a large number of experiments 
should be performed to determine the level of correlation between different end-point ob- 
servables. With a more detailed understanding of such correlations more constraints could 
be used than in the simple study described here, leading to further improvements in mass 
measurement precision. In addition it would be profitable to calculate analytical formulae 
for the end-point shapes similar to those for the global end-points described in Ref. M. 
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Appendix A: Three-body lepton-producing decay chains 

For completeness we list here some formulae for kinematic bounds equivalent to those listed 
in Section ||, for three-body lepton-producing decay chains incorporating single two-body 
quark-producing decays of the form: 

QL -> X° 2 Q Xihhq- (Al) 

The global end-points from such decay chains in NUHM models were considered in Ref. . 

The rriqu bound as a function of mu in this case is identical to that for the two-body 
case listed in Eqn. ( f4.3|) . The m q in^\ bound as a function of mu equivalent to Eqns. ( fP| ) 
and (fOj) can also be obtained from Eqn. ( |4.3| ) by recognising that m q in^\ is maximised 
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when the Xii q and one of the leptons are all co- linear in the rest frame of the x%- I n this 
case one of the two lepton+quark invariant masses is zero and hence 



K°(hO d(m " } ) 2 = Kr d(m " } ) 2 - ml (A2) 

_ ("4 - mio) 



2ml 

x 2 



(™|o - m|o + m\ ± ^A(m^o,m^o,m z/ )) . (A3) 



This is maximised when <m\\ = 0, which gives the position of the global end-point reported 
in Ref. 0: 

/ max \^ _ «f *2 y X 2 Xi> 

X2 



while if mjij = m™ ax = m^o — m^o then 



Xa ""Xi 

( max(m„)\2 _ V 1l x u 2 n X 2 x 2 Xi7 ,^ . 

I q/(hi) ) m 2 ■ v ^ 

X l 2 



The m g ;(io) bound equivalent to that derived from Eqn. ( |4.6f ) can be obtained from the 
three-body analogue of Eqn. ( |3.12 ) by conserving energy and momentum in the qll plane. 



This equation is more complicated than Eqn. (|3.12|) but setting m,ij = m q i 2 it simplifies 
considerably to give: 



b d( ) 2 (■^ / QL ^'~0 J / } j \ 

K^o" ) = 4m 2 X2 [ m l° ~ m \\ + m U ± V A (™*° » m X? » m «) j ( A6 ) 



fb°und(m a )\2 
l m g«(hi) J 



(A7) 



Consequently the m ? ^hi) bound is \/2 times as large as the m q i^ bound, as was noted for 
the global maximum in Ref. plf . 
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